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Abstract. We construct default-free interest rate models in the spirit of the 
well-known Markov funcional models: our focus is analytic tractability of the 
models and generality of the approach. We work in the setting of state price 
densities and construct models by means of the so called propagation prop- 
erty. The propagation property can be found implicitly in all of the popular 
state price density approaches, in particular heat kernels share the propaga- 
tion property (wherefrom we deduced the name of the approach). As a related 
matter, an interesting property of heat kernels is presented, too. 
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1. Introduction 

The heat kernel approach (HKA for short), which will be presented and discussed 
in this paper, is a systematic method to construct state price densities which are 
analytically tractable. To be precise on analytic tractability, we mean with this 
notion that bond prices can be calculated explicitly, and that caps, swaptions or 
other derivatives on bonds can be calculated up to one integration with respect to 
the law of the underlying factor Markov process. Therefore such models can be 
easily calibrated to market data and are useful for pricing and hedging purposes, 
but also for purposes of risk management. 

The original motivation of introducing the HKA was in modelling of int erest 



rates with jumps. In the HJM framework (jHeath. Jarrow and Morton! . 1 19921 ) . the 
drift condition b ecomes q uite complicated (see H. Sh irakawa's pioneering work 
( Shirakawa , 1991 ) , see also ( Akahori and Tsuchiva , 20061 ) and references therein) if 



taking jumps into account, while in the spot rate approach, one will find it hard 
to obtain explic it expressions of the bond pri ces (the affine class is almost the only 
exception). In ( Akahori and TsuchivaL [2006). the state price density approach^ is 



applied and by means of transition probability densities of some Levy processes ex- 
plicit expressions of the zero coupon bond prices are obtained. The HKA is actually 
a full extension of the method, and thanks to the generalization its applications are 
now not limited to jump term structure models as we will show in the sequel. 

Before the presentation of the theorem, we will give a brief survey of t he state 
price density approaches such as the potential approach by L.C.G. Rogers (Rogers], 
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1 The state price density, which is sometimes called the pricing kernel, or the state price defiater, 
is, roughly speaking, a Radon Nikodym derivative multiplied by a stochastic discount factor. A 
brief survey of the state price density and its relation to the interest rate modelling will be given 
in section UTT] 
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1997 ) or Hughston's approach ( Flesaker and Hughston , 1992f ) ( Lane P. Hughston 



20051 ), etc, in Section l2~Tl and Section l2~2l Our models are within t he class of Markov 



funct ional models proposed by P. Hunt, J. Kennedy and A. Pelsser (jHunt. Kennedy and Pelsser 



2000l ), and therefore compatible with their practical implementations. One of the 



contributions of the HKA to the literature could be to give a systematic way to 
produce Markov functional models. 

As a whole, this paper is meant to be an introduction to the HKA, together 
with some examples. Topics from practical viewpoint like fitting to the real market, 
tractable calibration, or econometric empirical studies, are not treated in this paper. 
Nonetheless we note that HKA is an alternative approach to practical as well as 
theoretical problems in interest rate modeling. 

2. State price density approaches 

2.1. State price density approach to the interest rate modeling. We will 

start from a brief survey of state price densities. By a market, we mean a family 
of price-dividend pairs (S l ,D l ), i £ J, which are adapted processes defined on a 
filtered probability space (O, T, P, {^t})- A strictly positive process (7r t ) t>0 is a 
state price density with respect to the market if for any i £ I, it holds that 

(2.1) Si = ^E[r \ s dD\\T t ], 

Jt+ 

or for any T > t, 

(2.2) S i t =Ttt 1 E[Tt T Si r + / w a dD\\F t ]. 

Jt+ 

In other words, T^T/^t gives a (random) discount factor of a cash flow 1 at time T 
multiplied by a Radon Nikodym derivative. 

If we denote by Pj the market value at time t of zero-coupon bond with maturity 
T, then the formula (12. 2| gives 



(2.3) Pj =^ x E\* T \T t \. 

From a perspective of modeling term structure of interest rates, the formula (|2.3j) 
says that, given a filtration, each strictly positive process 7r generates an arbitrage- 
free interest rate model. On the basis of this observation, we can construct arbitrage- 
free interest rate models. Notice that we do not assume TTt to be a submartingale, 
i.e. in economic terms we do not assume positive short rates. 

2.2. The Flesaker-Hughston model, the Potential Approach and an ap- 
proach by Wien er chaos. The ra tional log-normal model by Flesaker and Hugh- 
ston ( Flesaker and Hughston! . Il992h was a first successful interest rate model de- 



rived from the state price density approach. They put 

where A and B are deterministic decreasing process, W is a one dimensional stan- 
dard Brownian motion and a is a constant. The model is an analytically tractable, 
leads to positive-short rates, and is fitting any initial yield curve. Furthermore 
closed-form solutions for both caps and swaptions are available. Several extensions 
of this approach are known. 
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In (Roger sL ll99Tt> . L.C.G. Ro gers introduced a more systematic method to con- 



struct positive-rate models which he called the potential approach. Actually he 
introduced two approaches; in the first one, 

TTi = e- at R a g(X t )/R a g(X ), 

where R a is the resolvent operator of a Markov process on a state space S, g is a 
positive function on S, and a is a positive constant. In the second one, g is specified 
as a linear combination of the eigenfunctions of the generator of X. Note that when 
X is a Brownian motion, e for any 9 £ R is an eigenfunction of its generator, and 
the m odel gives another perspe ctive to the rational lognormal models above. 



In (jLane P. Hughstonl , 120051 ) L. Hughston and A. Rafailidis proposed yet an- 
other framework for interest rate modelling based on the state price density ap- 
proach, which they call a chaotic approach to interest rate modelling. The Wiener 
chaos expansion technique is then used to formulate a systematic analysis of the 
structure a nd classification of interest rate models. Actually, M.R. Grasselli and 



T.R. Hurd (Grasselli and Hurd. , 2003) revisit the Cox-Ingersoll-Ross model of in- 



terest rates in the view of the chaotic representation and they obtain a simple ex- 
pression for the fundamental random variable as 7r t = -E^X^ — (E[XaD | Tt]) 2 \ Tt] 
holds. In line with these , a stable property of the second-order chaos is shown in 
(jAkahori and Haral . [20061 ) . 



3. The Heat Kernel Approach 

In this section we present the simple concepts of the heat kernel approach (HKA) 
and several classes of possibly interesting examples. The guiding philosophy of HKA 
is the following: if one can easily calculate expectations of the form E(f(X t )) for 
some factor Markov process X, and if one knows explicitly one additional function p 
depending on time and the state variables, then one can construct explicit formulas 
for bond prices, and analytically tractable formulas for caps, swaptions, etc (i.e., 
those formulas are evaluated by one numerical integration with respect to the law 
of the underlying Markov process). In other words, any explicit solution of a well- 
understood Kolmogorov equation with respect to some Markov process leads to a 
new class of interest rate models. 

3.1. Basic concept. Wc consider a general Markov process (X x ) t>0 x£S on a pol- 
ish state space S and a probability space (fi,^ 7 , P) with filtration (^Ft) t >o- 

Definition 3.1. Letp(t, x) defined on R>o x S be a strictly positive function. It is 
said to satisfy the propagation property with respect to X if 

(3.1) E(p(t,X*))=p(t + s,x) 

for all t > 0, s > and x £ S. We assume furthermore that p(t,X t ) £ L 1 (f2) for 
t > 0. 

Let p satisfy the propagation property with respect to the Markov process X 
and p(0,x) £ L 1 (R). We define then a state price density Wt through 

(3.2) Tr t :=p(t,Xf). 



We clearly have the following basic assertion: 
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Theorem 3.2. Bond prices are given through 

pT ^p{2T-t,Xf) 



p{t,Xf) ' 

for < t < T and x G S . 

Proof. The bond prices are calculated with respect to the state price density 
through 

p t = ±Efa\r t ). 

Actually, by the Markov property and by the propagation property we can calculate 
the right hand side, 

E(ir T \Ft) = E(p(T t XZ)\F t ) 

= E(p{T,X%._ t ))\ y=x . 
= p(2T-t,Xf) 

for all x € S and < t < T. □ 

Remark 3.3. Notice that this remarkable simple formula is given with respect to 
the historical /physical measure. 

A derivative with payoff G on Pj for some < t < T has at time t = the price 

-E[p(t,X t )G( )j, 

which is analytically tractable if one knows how to calculate expectations of the 
form£(/(Xf). 

Remark 3.4. Notice that the initial term structure is given in the previous setting 
through 

P( T) = P(2T ' X) 
l ' j P(0,x) 

forT> 0. 

3.2. Generic example. Let {Xf ) t>0 IgS with a polish state space S on a proba- 
bility space {fljJ 7 , P) with filtration (J 7 t) t>0 . Take a measurable, bounded h : S — > 
R>o, then 

(3.3) p(t,x):=E(h(X?)) 
satisfies the propagation property (|3.ip . 



This example demonstrates that the previous theory covers the results of (lAkahori and Tsuchiva , 



20061 ) . Indeed, consider a d-dimensional Levy process (X t ) t>0 starting at in its 



natural filtration. We assume that it has a density y h- > p(t, y) with respect to 
Lebesgue measure. Then Xf := x — X t defines a Markov process on M. d and we 
have that 

(3.4) E(p(t,x-X s ))=p(t + s,x) 

for s, t > and i G M d . 

Notice that this example can be associated to the previous Example 13.21 if one 
allows h = 5q in the sense of distributions. 
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Figure 1 . Simulated Yields in the generic affine class driven by CIR 



3.3. A generic affine example. Consider an affine Markov process X on R™ x 
R m and m, ...,fi r e R n+m such that E(exp((jM,X t ))) exist. Then we define ~ 

r 

h(x) = J^ai exp((fij,x))) 

z=l 

for some numbers Oj > 0. In the spirit of the previous generic example we obtain 

r 

P(t, x) = ^2 <kE(exp{{lH, X t ))) 
i=i 

for t > and each initial value x £ K" x R m . Notice that the "Laplace transforms" 
E(exp((ni, Xt))) are easily calculated for affine processes, which yields a concrete 
representation of the prices in the physical measure. These models are completely 
new interest models constructed from affine models. For instance, the derived 
short rate models have interesting (econometric) features. An example driven by 
the affine Cox- Ingersoll- Ross process is visualized in Figure 1. 

3.4. Pricing formula of Swaptions. The pay-off (at its maturity T a ) of a (unit 
of) swaption is 



(3.5) (S a , (T a ) - K)+ nP(T a ,T t ), 

where {T Q , T Q+ i, . . . , Tp\ with r, = Tj — Tj_i is the tenor structure of the swap 
contract, and 

\-P{T a ,Tp) 

is the swap rate (see e.g. ( Brigo and Mercuric! l2006l )*Chapter 1, section 6). Note 
that (13.51) can be rewritten as 



(3.6) {l-P(T a ,Tf})-K J2 nP(T a ,T l )}+. 

i— a+1 
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Therefore in the state price density approach, the quantity 

1 
SWt :=-E(TT Ta {l-P(T a ,T )-K V TiP(T a ,Ti)} + \r t ) 

gives a fair price of the swaption. Substituting 

P(T a ,T i ) = E(-K Ti \F T J/ir Ta , 

we get 

1 13 
sw t = — E({ir Ta - E^t^TtJ -K nE(ir Tl \FT a )} + \Ft). 

i=a + l 

By HKA, namely 

7T t =p(t,X t ), 

where p has the propagation property w.r.t. a Markov process X. We have 

1 



p(t,X t ) 

(3. 7 ) x £ Xt [{p(T„,lT a - f ) -P(2T -T Q ,X Ta _ t )- 

-if ^ 7-^(21; -T a ,Xr -t)} + ], 
which is again up to integration with respect to the law of X - explicit. 



3.5. Relations to short rate models. Assume that p(t, x) is differentiate with 

t 



respect to t. Then the short rate Rt — —dr log Pf\t=T of the previous interest rate 



model is given through 

^ = -p^xf)Wt^ x ^ 

Since p > the sign of the short rate is determined through 

d_ 



--p(t,X*). 



Notice that we do not necessarily have a short rate process. 

3.6. Eigenfunction models. Take a finite dimensional manifold M and take X 
a Markov process with values in M. A non- vanishing function g : M — > R with 

(3.8) E{g{Xf)) = exp^t)g(x) 

for some fi < is called eigenfunction of X. Let : M — > M, i = 1, • • • , JV be 
eigenfuntions with eigenvalues /i, < 0, and let Ai : M>o — > K>o be decreasing 
functions. Set 

N 

v(t,x) := Ai(t)5i(a;), 

i=l 

then we notice that v(t, X t ) is a supermartingale, too. If v is strictly positive 
and N > 1, then we can construct an arbitrage- free positive interest rate model by 
setting ir t := A t ) as a state price density. We call such models the eigenfunction 
models. 



Remark 3.5. The class includes the generic approach II in the potential approach 
by Rogers tRoaers\ . \l99't) , and is included in the Flesaker-Hughston's general frame- 
work since e^^ t g{X t ) becomes a martingale (see section \2. l 2\) . 

3.7. Swaption Formula via eigenfunction models. The formula (|3.7[) becomes 
extremely simple when we use the eigenfunction model. 

(3.9) n =p(t,X?) = l + e>*g(X?) 

where g satisfies (I3.8|) . By applying <\'S.9\i on (|3.7p , we obtain the following result. 
Theorem 3.6. The pricing formula of the swaption is 

where 

A = e »T a _ e ri2T f ,-T a )_ K J- r . e ^m+^-T a)) 

i— q+1 

and 

B = K(Tp — T a ). 

Example 3.7 (Brownian motion with drift). Let X be a Brownian motion with 
drift: whose generator is A/2 — k ■ V where n € M. d . For g(x) = e^ c ' x ^ with c £ R d ; 
we have 

(A/2-K-V)g=(\c\ 2 /2-(K,c))g. 

Therefore, if c is such that \c\ 2 /2 — (k,c) < (this is always possible, take c = k 
for example), then the state price density 

p(t,X t ) = l + e"g(X t ), 

where v :— |c| 2 /2 — (k,c), gives a positive rate model. In this case 

where 

BS(a, b, c, d) = [ (ae cx+d - b) + e~^ dx. 

V 27T Jffi 

Notice that the Black-Scholes formula for the price at time t of the call option with 
the strike K and the maturity T is given by 

e -r(T-t) BS ( St ^ K ^ a ^r~t, (r - °— )(T - t)). 

Remark 3.8. This example corresponds to the Flesaker-Hughston's rational log- 
normal model \Flesaker and HuahstoA. \l99& ). 

Example 3.9 (squared OU process). Let X be an Ornstein-Uhlenbeck process on 
W 1 , whose generator A is given by 

A := -fix\7 + 

where p < 0. Setting g(x) = e^ x \ , we have 

Ag{x) = pdg(x) 
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Therefore we have that g is an eigenfunction with the eigenvalue fj,d. Setting d = 1 
and pit, x) = 1 + g[x), we obtain 

(3.10) SWt = Ae^-t)g{X?){ <$>{d+) - B{<S>(d+) - *(^)} 



1 + e^g{Xf) 



where <&(x) — f e ,„ dx. 



± V » 



± J I log | - Xf e M(ra-t) 



(3.11) df = 

V(2M)- 1 (e MTo_t) - 1) 

and 



iJllogf-Xfe-^-*) 



(3.12) d± - . 

N /(2/x)- 1 (l-e- 2 ^ T --*)) 

For i/ie proof see the Appendix\5\ 

4. HKA and Positive Rate Models 

As we have pointed out the propagation property of p does not ensure that 
the process p(t, X t ) is a supermartingalc. Thus, the interest rate model with ir t = 
p(t, X t ) does not exclude the possibility of negative rates. In this section we provide 
models based on HKA and leading to positive short rates. 

4.1. Weighted heat kernel approach. Let p be a positive function with propa- 
gation property with respect to a Markov process X, and a weight / : K>o x K>o — ¥ 
M>o such that fit, u — s) < f{t — s, u) for arbitrary t,u G M>o and s <t Au. Set 

/>oo 

q(t,x) = / p(s,x)f(t,s)ds, 
Jo 

Then we have 

Proposition 4.1. The process q(t,X t ) is a supermartingale. 
Proof. For t > s and x € S, we have 

/>oo 

E{q{t,Xf)\r.) = / £(p( u ,A7)|.F s )/(^)d w 



p(u + 1 — s, Xf)f(t, u) du 
p{v,X*)f(t,v-t + s)dv 

3 

Here we have used the propagation property of p. Since 

i poo 

p(v, Xf)f(t, v-t + s)dv< / p(v, Xf)f{t - (t - s), v) dv 



t-s 



< / p(v,X*)f(s,v)dv = q(s,Xf), 
Jo 

we have the desired result. □ 



Thus interest rate models with Tr t = q(t,X t ) ha ve posit i ve sho rt rates. The class 
is an extension of the potential models defined in (jRoger l l2006lh since we can take 
p as in (O and f{t, s) = e - a( - t+s K 



4.2. Killed Heat kernel approach. Let V be a non-negative, measurable func- 
tion denned on S. Put 

ft 



q(t,x)=E(exp(- f V(X*)ds)), 
Jo 



and 

(41) p r =q (2T-t,Xn 

Then the bond market spanned by P T is arbitrage-free with respect to the state 
price density 



ir t = q(t,X?)eM- I V(X*)ds). 
Jo 



Indeed, this corresponds to an additional factor = y + J Q V(Xf ) ds for y > 
added to the Markov process X. Since (X X ,Y V ) is again a Markov process on 
S x M>o, we notice that 

p(t, x, y) := E{e- Y ?) = e -"f?(exp (- / V(X*) d*)) = e~ v q{t, x) 

Jo 

defined on R>o x S x R>o satisfies the propagation property. Therefore the previ- 
ously described bond prices are arbitrage-free with respect to the stated state price 
densities for initial values (x, 0). Since q is decreasing in t the short rates are again 
positive. 

4.3. Trace approach. We give a third, more involved method to obtain positive 
rate models. Though it is only valid for symmetric Levy processes at this stage, 
again the class is completely new. Let (X t ) t>0 be a d-dimensional Levy process, 

which is symmetric in the sense that X x = -X 1 . We define X? = x + X t . Note 
that in this case its Levy symbol is positive and symmetric, i.e. 

for some 

ip : R d -> M> , 

which is symmetric in the sense that VKO = VK - 0- 
Let /j, be a measure on R d such that 

/x(-A) = fj,(A) 

for any Borel set A. Put 

u(t,x)= [ e^^^e-^^drj), 



assuming it is finite (except at t = 0) for each x € R d . 

Lemma 4.2. The function u is non-negative and has the propagation property with 
respect to {X x }. 

Proof. First, notice that by the symmetry of \x we have 

(4.2) u(t,x)= [ cos(2Tr(n,x))e- t ' l ' {v) ii(dn) e R, 



and u(t, x) — u(t, —x) for any x G Mr. Then by Bochner's theorem we know that 
u(t, •) is non-negative definite and we see that u(t, x) > for any x £ Mr. 
Since < 1, E[u(s, Xf)] < 1 for any s, t > and x e M d , we have 

E{u{s,X?)) = E (j e^to' x Ve-*Wp(d<n)j 

□ 

Remark 4.3. J/ i/ie Fourier transform 

f(x) := / e 2 "<^>/x(dO 



exists and is bounded in x, then by the symmetry of p, it holds that fix) > ; and 
u has the following expression: 

u(t,x)=E(f(X?)). 

On the other hand, let p be the Lebesgue measure of M d . If the density pt (x) of X® 
exists in L 1 n L? , we have 

u(t, x) = p t {x). 

That is, u is the heat kernel. 

We are now able to prove the following theorem on supermartingales constructed 
from heat kernels: 

Theorem 4.4. Let c be a positive constant greater than 2 and let A be any positive 
constant. Set n t :— u(A + t,Xf) + cu(X +£, 0). Then 7r is a supermartingale with 
respect to the natural filtration of X . 

Proof. By the expression (|4.2p , we have 

u(A + 1, Xf) - u(X + 2T-t, X?) 

(43) =/ cos2 7 r{Xf,0(e-( t+A )^«)-e-( A+2T - 4 )^(f)) M (dO 



Since iji > and A is positive, 
RHS of i 

> f ( e -(W-A)*«) _ e -( A + 2T -*)W«))^) 
= + A, 0) - u(A + 2T - 1, 0)). 

Thus we obtain 

u(t + A,X t x ) + it(t + A,0) 

> it(A + 2T-t, Xf) + u(A + 2T-t, 0), 
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and then 

u(t + X,Xf) + cu(t + X,0) 

(4.4) 

V 1 >u(X + 2T-t,X?) + u{X + 2T-t,Q) + (c-l)u{t + X,Q). 

On the other hand, by an elementary inequality we have 
u(X + 2T - t, 0) + u(t + A, 0) 

(e -(A + 2T- tWe) + ^(t+A)^))^) 

i 

> 2 / e-^ A + A + 27 >«V(^), 

and since A is positive, we obtain 

u(X + 2T - t,Q) + u(t + A,0) 

(4 ' 5) > 2 / e-( T+A ^«V(^) - MT + A, 0). 



Combining (|4.4|) and (14. 5|) , we have 

u(t + A,Xf) + cu(t + A,0)- 
(4.6) -u(A + 2T-i,A: t a; ) -cu{T + A, 0) 

> (c - 2)(u(i + A, 0) - (c - 2)u(T + A, 0)). 

The right-hand-side of (|4.6p is non-negative since c > 2 and 

u(t,0)= / e-'«SK) 
is decreasing in i. Finally, since we have 

^M-? 7 *] = u(A + 2T-t, X'f) + cu{T + A, 0), 
we now see that wt is a supermartingale. □ 

4.4. Examples. By the above theorem, we can now construct a positive interest 
rate model by setting u(X + 2T — t, Xf) + cu(T + A, 0) to be a state price density, 
or equivalently 

pT _ u(X + 2T-t, Xf) + cu{T + A, 0) 

* ~~ u(t + X, X?) + cu(t + X, 0) 

Here we give some explicit examples. In all examples, c is a constant greater than 
2 and A some positive constant. 

Example 4.5 (Heat Kernel with Trace). A simple example is obtained by letting 

1 

u(t, x) = e 2* 

V2^i 

where the bond prices are expressed as 




Figure 2 . Simulated Yields for the trace class model driven by a 
Quadratic Gaussian 



Example 4.6 (Quadratic Gaussian with Trace). A similar model can be obtained 
by setting, for a > 0, 



u(t, x) 



1 



_g 2(ta^+l) 



= E 



e 2 dy 



A visualization of this model is shown in Figure 2. 

Example 4.7 (Cauchy with Trace). Let Z x be a Cauchy process in M d starting 
from x. The explicit form of the transition density of Z° is known to be 

T((d+l)/2)9t 



(4.7) 



P(t,x) 



(irKet^ + lx-t^})^ 1 )/^ 

where 9 > and 7 G Mr (see e.g. Note that Cauchy processes are 

strictly stable (or self-similar) process with parameter 1; namely, 

Z a t = aZ t , a> 0,t > 0. 
This property is seen from their Levy symbol: actually we have 

E[e i{ ^ Zt} ] = e -*( 9 lfH<«>7». 

If we set 7 = 0, the density turns symmetric and therefore, setting 7r t = p(X + 
t,Zf) + cp(A + t, 0), we obtain an explicit positive term structure model. We give 
the explicit form of the bond prices in the case of d = 1: 

8 2 (\ + t) 2 + \Z?\ 2 



P 



X + t 

A + 2T- t 9 2 (\ + 2T- t) 2 + \Zf\ 2 

9 2 (l + c)(\ + 2T-t) 2 + c\Zf\ 2 



9 2 {\ + c){\ + t) 2 +c\Z'f\ 2 ' 
The model is a modification of the Cauchy TSMs given in {Ak ahori and Tsuchivd . 



2006 
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Let X t be Levy process and T t be a subordinator independent of X t . Recall that 
the process Z t := Xt ± is usually called the subordination of X by the subordinator 
T . Among the subordinations of the Wiener process, two classes are often used in 
the financial context; one is the Variance Gamma processes and the ot her is Norma l 
Inver se Gaussian processes, since they are analytically tractable (see ([Rama Contl . 
120041 ) and references therein). The subordinator of a Variance Gamma process 
(resp. Normal Inverse Gaussian process) is a Gamma process (resp. Inverse Gauss- 
ian process). If we let the drift of the Wiener process be zero, then the densities of 
the subordinated processes satisfies the condition for Theorem 14.41 

Example 4.8 (Variance Gamma TSMs with Trace). The Gamma subordinator T G 
is a Levy process with 

P(T t G G dx)/dx = ^ x *-l e -^l {xm , 

where r\ and 7 are positive constants. Then the heat kernel of the subordinated 
Wiener process W t g is 

(4.8) u(t,x):= e -^J— s ^e^ s d S . 

Jo V^ns T(r]t) 

From this expression we have 

The heat kernel {4.8\ ) can be expressed in terms of the modified Bessel functions 
K p , p G M, which has an integral representation 

(4.9) K p (x) = ±(^y J°° e-t-^t-v-Ut, x>0. 
Actually we have 

U (t,x) = J^mr 1 (iy^ +i \x\-^K_ 11t+h (\x\^¥f 



Example 4.9 (Normal Inverse Gaussian TSMs with Trace). The heat kernel of 
the inverse Gaussian subordinator T/ is 

X ' 

where r\ and 7 are positive constants. The heat kernel of the subordinated Wiener 
process W T i is 

u(t 7 x) := / e~^pt(s)ds. 
Jo v 27TS 

Applying J^.g[ ), we obtain 

u(t,x) = ^-^e^^Ki (^7(2^ + 4^2)) . 
In particular, we have 

u(t,0) = ^^Le^'^K, (2rfofw) . 
2?r 2/07 
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5. Appendix 

Lemma 5.1. For m,a,b € M, v > 0, and fi < 0, it holds 



(5.1) 



Proof. 



2 1 (x-m) 2 

e Ma: e 2 « da; 

,, ., i, \/2ttv 

„ j-i/ .,.a-m. 

_ e l-2(.« $( _) - $( _). 



w/iere 

(5.2) 77i := and w := 



1 — 1 — 2/ii> 



2 1 (X-TT1) 2 

e nx e 2 « dec 

a<a;<b \j2tTV 



V k-™ 2 f 1 O-m) 2 

(5.3) =-i/-e 1 - 2 M« / e 2 * da; 

V « ./a<x<& V^TTV 

V _ £i m 2 _ / ,5— 777 , ,a™ 777 N \ 



A proof of ([3^0]) with (pLl!]) and (|37L2~1) . Notice that 

(5.4) x T a -t\ Xt ~ N{m u v t ) 
where 

(5.5) m = X t e-^ T "-V and u = — (1 - e" 2 ^ 
Put 



a = -\/^ log f and 6 = \/^ log T 

By Lemma l5.il we have 

f (Ag(y) - B)+P(X^_ t e dy \ x = X t ) 

( 5 6) 

ffj „ m 2 ( b — m. ,.a — m. \ / ' 6 
= ii -e'-v $( — r=^) - $( — — ). - B\ $(— 



By (jOj) and ([531) 



e 



2// 



and 

With these expression we have (|3.10|) . 
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